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ABSTRACT 

We report interferometric mapping of the bipolar pre-planetary nebula IRAS 08005-2356 (108005) 
with an angular-resolution of ~l”-5”, using the Submillimeter Array (SMA), in the ^^CO J=2-l, 

3-2, ^^CO J=2-l and SiO J=5-4 (v=0) lines. Single-dish observations, using the SMT 10-m, were 
made in these lines as well as in the CO J=4-3 and SiO J-6-5 (v=0) lines. The lines profiles are very 
broad, showing the presence of a massive (> 0.1 A/q), extreme high-velocity outflow {V ~ 200 km s“^) 
directed along the nebular symmetry axis derived from the HST imaging of this object. The outflow’s 
scalar momentum far exceeds that available from radiation pressure of the central post-AGB star, and 
it may be launched from an accretion disk around a main-sequence companion. We provide indirect 
evidence for such a disk from its previously published, broad Ha emission profile, which we propose 
results from Ly/3 emission generated in the disk followed by Raman-scattering in the innermost regions 
of a fast, neutral wind. 

Subject headings: circumstellar matter - planetary nebulae: individual (IRAS 08005-2356) - accretion 
disks - stars: AGB and post-AGB - stars: mass loss - stars: winds, outflows 


1. INTRODUCTION 

Following the ejection of half or more of their mass 
via isotropic, slowy expanding winds, AGB stars evolve 
into planetary nebulae (PNe), which, surprisingly show 
a diverse range of aspherical (e.g., bipolar and multi¬ 
polar) morphologies (e.g., Sahai & Trauger 1998, Sahai, 
Morris & Villar 2011a). Studies of pre-Planetary nebu¬ 
lae (PPNs), objects in transition between the AGB and 
planetary nebula (PN) evolutionary phases, are critical 
for characterising the physical processes responsible for 
this dramatic transformation. Sahai & Trauger (1998) 
proposed that fast collimated outflows or jets, operat¬ 
ing during the PPN and/or very late AGB phase, are 
the primary agents for the dramatic change in the mass- 
loss geometry and dynamics during the AGB-to-PN evo¬ 
lutionary phase. However, the physical mechanism for 
producing these fast outflows remains a mystery. High- 
angular-resolution interferometric (sub)millimeter-wave 
observations are the best way to quantitatively probe 
the fast outflow’s dynamics and energetics - crucial in¬ 
formation for theoretical models (e.g., Akashi & Soker 
2013), and detailed numerical hydrodynamical simula¬ 
tions (e.g., Lee & Sahai 2003, Balick et al. 2013), for PN 
shaping. 

Such observations have resulted in the discovery of 
a handful of “extreme-outflow” PPNs - objects in 
which the molecular outflows reach speeds in excess 
of ~100kms“^, e.g.. Boomerang Nebula (Sahai et 
al. 2013), IRAS 22036-f5306 (Sahai et al. 2006), 
IRAS 19374-1-2359 (Sanchez Gontreras et al. 2013/ 
IRAS 16342-3814 (Imai et al. 2012) and HD 101584 
(Olofsson et al. 2015). Detailed studies of such ex¬ 
treme objects are likely to provide the best motivation 
for, and most stringent tests of, theoretical models to 
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explain their origin (e.g., Blackman & Lucchini 2014). 

In this Letter, we report (sub)millimeter-line obser¬ 
vations of the PPN IRAS 08005-2356 (108005), which 
clearly reveal it to be an extreme-outflow PPN. Early CO 
J=2-l observations by Hu et al. (1994) resulted in the 
marginal detection of a weak, broad line. 108005’s F51e 
central star V510Pup (Slijkhuis et al. 1991) may have 
made a recent transition from ejecting oxygen-rich mate¬ 
rial to carbon-rich material (Bakker et al. 1997). Its mor¬ 
phological classification is Bo*(0.55) (Sahai et al. 2007), 
i.e., it has a bipolar morphology (resolved via HST imag¬ 
ing, Ueta et al. 2000) with lobes open at their ends and a 
central star seen at 0.55/rm. Its optical spectrum reveals 
the presence of a prominent Ha emission line with very 
broad wings (FWZI~ 2400 km s“^) and a P-Cygni type 
blue-shifted absorption feature (Sanchez Contreras et al. 
2008: SCetal08; also Slijkhuis et al. 1991, Klochkova & 
Chentsov 2004). Its estimated distance ranges between 
2.85 kpc (Oppenheimer et al. 2005: OBS05) and 3-4 kpc 
(Klochkova & Chentsov 2004); we adopt a value of 3 kpc. 

2. OBSERVATIONS 

The 1.3 mm and 0.87 mm interferometric observations 
were obtained with the Submillimeter Array (SMi^ 
at Mauna Kea, Hawaii. Bandpass calibration was per¬ 
formed using observations of 3C279. At 1.3 (0.87) mm, 
complex gain calibration was obtained from observations 
of the quasars 0750-1-125 and 0730-116 (0747-331 and 
0826-225), and flux calibration was obtained from obser¬ 
vations of Callisto (EuropaL Additional observing pa¬ 
rameters are listed in TableflJ 
SMA data were calibrated using the MIR-IDL pack- 
age[J and images were made using the Miriad software. 
Data cubes were obtained with a velocity resolution 

® The Submillimeter Array is a joint project between the Smith¬ 
sonian Astrophysical Observatory and the Academia Sinica Insti¬ 
tute of Astronomy and Astrophysics and is funded by the Smith¬ 
sonian Institution and the Academia Sinica 
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smoothed to 10 km s“^ per channel (to increase the S/N 
in each channel). Natural weighting was used to produce 
all images. 

Single-dish observations of the J=2-l, 3-2 and 

4-3, J=2-l, and SiO (v=0) J=5-4 and 6-5 line 

emission were obtained at SMT during 2014, Nov/Dec 
and 2015 Jan. Telescope pointing was frequently checked 
on VY CMa, and is estimated to be better than a small 
fraction of the beam. The weather was generally good, 
with system temparatures in the range Tsys 200 — 
250 K at 1.3 mm and Tgyg ^ 1300 K at 0.7 and 0.8 mm. 
Linear baselines were subtracted from the spectra shown. 
We have assumed main-beam efficiency factors of 0.76, 
0.66, and 0.62 for the 230.6, 345.8, and 461.0 GHz obser¬ 
vations (Edwards et al. 2014). 

The optical images were taken with the NASA/ESA 
Hubble Space Telescope in GO programs ID 6364 (PI: 
M. Bobrowsky) and 6366 (PI: S. Trammell), using the 
PC on WFPC2, and extracted from the Hubble Legacy 
Archiv^il 

3. RESULTS 
3.1. Optical Imaging 

108005 shows an hourglass morphology in the HST 
F439W, F555W, and F675W images (Fig.|H - the lobes, 
separated by a narrow dark lane, flare out from the cen¬ 
tral star’s location, and attain a cylindrical shape. The 
SE lobe is significantly brighter than the NW one. A 
narrow, slightly-curved feature of linear extent 0738, is 
seen extending from the central star in the SE lobe (see 
inset). The feature is oriented at a PA = 120°, slightly 
different from the PA of the nebula symmetry axis, 143°. 
The total linear extent of the nebula along its long axis, 
as seen in the F555W image, is ^ 276 (1.2 x 10^^ cm). 

3.2. Millimeter-Wave Observations 

We detected the CO J=2-l, 3-2 and SiO J=5-4 lines 
with the SMA. The CO J=3-2 and SiO J=5-4 have sig¬ 
nificantly lower signal-to-noise, and the ^^CO J=2-l is 
tentatively detected. The (spatially) integrated SMA 
CO J=2-l line profile (using a 15" diameter circular 
aperture) shows a very broad profile (Fig.I^ covering 
a total velocity extent of about 350 km s“ , i.e., from 
-150 ^ Visr (km s-i) ^ 200. The SMT CO J=2-l line 
profile (Fig.[^ is similar to the SMA one (a narrow emis¬ 
sion feature at Visr = 36 km s“^due to the presence of 
an unrelated line-of-sight interstellar cloud, has been re¬ 
moved from the profile by interpolation). A Gaussian fit 
to the profile gives a central velocity Visr ~ 25 km s“^, 
and a FWHM of 196 km s“^. 

We have divided the velocity range spanned by the 
CO J=2-l line into three intervals on each side of 
the central velocity - the extreme high-velocity or EV 
component (blue: —150 < Visr (km s“^) < —85, red: 
130 < Visr (km s“^) < 200), a medium-velocity or MV 
component (blue: —85 < V)sr (km s“^) < —30, red: 
70 < Visr (km s“^) < 130), and a low-velocity LV com¬ 
ponent (blue: —30 < V;sr(kms“^) < 25, red: 25 < 
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Visr (km s“^) < 75). The SMA CO J=2-l emission in 
the blue and red-shifted parts of the EV, MV, and LV 
components show separations of about 271, 171, Oi'S, re¬ 
spectively, along the nebular symmetry axis (Fig.|^). 

These results imply the presence of a fast, bipolar out¬ 
flow with a (roughly) linear velocity-gradient, directed 
along the nebular symmetry axis, and a linear extent 
of ~ 271. Since AGB mass-loss is typically spheri¬ 
cal, with outflow velocities of 10-15 km s“^, we conclude 
that the outflow observed in CO J=2-l emission is a 
fast, collimated post-AGB outflow, typical of the PPN 
evolutionary phase. We made similar plots from the 
SMA data cubes for CO J=3-2 (Fig.|^) and SiO J=5- 
4 (not shown), dividing the emission velocity-range into 
two halves only (blue: —90 < V)sr (km s“^) < 20, red: 
30 < Visr (km s“^) < 130) due to the lower S/N in these 
lines, and found that these show a separation of 0785 
roughly along the nebular axis. This separation is equal 
to that derived from our CO J=2-l data for the same 
red and blue outflow-velocity ranges, implying that the 
SiO outflow’s linear extent is likely comparable to that 
of CO. 

We also detected the CO J=4-3, J=3-2 and ^^CO 
J=2-l as well as the SiO (v=0) J=5-4 and 6-5 lines 
with the SMT, but with lower S/N (Fig.[^. We have 
compared the total CO J=2-l line fluxes from our SMA 
and SMT data, as these have the highest S/N, and find 
that these are consistent within 15%, i.e., within typical 
calibration uncertainties) - we conclude that the SMA 
data do not suffer from any significant flux losses. 

Continuum images were obtained after removing spec¬ 
tral regions containing line emission, using the Miriad 
task uvlin. No continuum emission was detected from 
108005. Using line-free channels, we find a Itr noise of 
1.05 (1.2) mJy/beam in the USB (LSB) at 1.3 mm and 
11 (18) mJy/beam in the USB (LSB) at 0.87 mm. 

4. THE POST-AGB BIPOLAR OUTFLOW 
4.1. Outflow Properties 

We now determine the physical properties of the high- 
velocity bipolar outflow (e.g., scalar momentum, mass, 
and age) from an analysis of the molecular-line data. 
The observed ratios of the CO J=2-l, 3-2 and J=4-3 
fluxes (in Kkm s“^) are 1 : 1.27 : 1.2, which imply cor¬ 
responding non-beam diluted flux ratios of 1 : R 32/21 : 
R- 43/21 = 1 : 0.6 : 0.3, since the source is unresolved by 
the SMT beams of 32", 22", and 16" in these lines. 

We first make a simple emission model assuming a 
spherical outflow at a uniform excitation temperature 
equal to the kinetic temperature, ^in, and carry out 
a least-squares fit to the line fluxeR varying Tkin, the 
outer radius (Rout), and mass of the emitting region 
(Mg). We find Rout = 1-4" and Tki„ = 13.6K. How¬ 
ever, the derived Tkin is rather low compared to energy 
of the CO J=4 level (55.4 K), and modeling with the 
non-LTE RAD EX code (van der Taak et al. 2007) shows 
that the J=4-3 excitation temperature, Tex( 43 ), is lower 
than Tkin by about (20-25)% due to the relatively low 
CO line optical depths in our model (< 0.01) and the 
average density, Uav ^ 3.5 x lO'^cm”^, implied by the 
emitting region’s mass and volume. Thus, the resulting, 

® derived using the formulation in Olofsson et al. (1990) 
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corrected model -R 43/21 0 . 11 , too low to fit the data. 

Using RADEX, we find that the average density re¬ 
quired to bring Tex( 43 ) closer to Tkin in order to fit 
the observed source brightness temperature ratios, is 
Uav > 7.5 X lO'^cm”^; the required Tkin = 15.5K, giving 
a source size of 274, comparable to the CO source-size 
estimated from the SMA data. The model source size is 
similar to the size derived from the HST image, implying 
that the CO emission likely comes from the dense walls 
of the bipolar lobes seen in scattered light. The CO col¬ 
umn density is, Nco > 5 x 10^^ cm“^ (and the J=2-l, 
3-2, and 4-3 lines have optical depths of 0.48, 0.44, 0.20). 
Assuming that the average emitting column is equal to 
the source radius, we find a CO-to-H 2 abundance ratio of, 
fco = 1.3 X 10“^, in reasonable agreement with the value 
typically assumed for PPNs, 2 x 10“^ (e.g., Bujarrabal 
et al. 2001). Assuming a spherical emitting volume, the 
mass is. Mg = 0.076 M©. The ^^CO/^^CO abundance ra¬ 
tio is, f(i2C/i3C)=9.6, from fitting the i^CO/^^CO J=2- 
1 line flux ratio corrected for the different beam-dilutions 
( 8 . 1 ). 

The above values of Mg and f(^^C/^^C) are lower limits 
since this is our “minimum-mass” model - models with 
higher values of nav are allowed. However, the total mass 
increases more slowly than the average density, since in 
models with higher values of Uav, the emitting region’s 
size is smaller. For example, if we use the dust mass, 
Md=O.OO19L0 (Zl/2.85kpc)^ derived by OBS05 from a 
detailed 2D radiative transfer model of 108005’s SED, 
scale it to D=3 kpc, and adopt a typical gas-to-dust ra¬ 
tio for oxygen-rich AGB stars, Mg/Md = 200, we get 
Mg = 0.42 Mq. For this value of Mg, we need a CO 
model with nav ~ 8 x 10 ®cm“^, a factor ~ 10 higher 
than in our minimum-mass model; Tkin = 11.5 K, and 
^CO = 1-6 X 10^®cm“^. Since the CO J=2-l opti¬ 
cal depth is higher (t 2 i = 1.95), the abundance ratio 
is higher, . 

We calculate the scalar momentum using the formu¬ 
lation described in Bujarrabal et al. (2001). Using our 
minimum-mass model, we find Pgc 2.8 x 10 ^® gcms“^ 
for an inclination angle of the nebular axis to the sky- 
plane, i = 30° (OBS05). The kinetic energy in the 
outflow is Ekin ~ 2.6 X 10^^ erg. These values of 
Pgc and Ekin He near the upper end of the range for 
PPNs, gcms“^ and 10^®““^®erg (Bujarrabal et 

al. 2001). This outflow cannot be driven by radia¬ 
tion pressure because the lobes’ dynamical (expansion) 
time-scale, tdyn = 190 yr (from dividing the model CO 
shell-size by the CO J=2-l FWHM line-width) is much 
smaller than that required by radiation pressure to accel¬ 
erate the observed bipolar outflow to its current speed, 
trad = Pscl{Ljc) ~ 6.6 X 10^ yr, given 108005’s lumi¬ 
nosity of 6980 Lq at D=3kpc (using the value derived 
by OBS05, 6300(11/2.85kpc)® L©.) The mass-loss rate 
in the outflow is > 5.8 x 10“"^ M© yr“^. 

The observed SMT SiO J=6-5 to 5-4 line flux ratio 
is 0.94 (with about ±15% uncertainty), implying an in¬ 
trinsic flux ratio (i.e., corrected for beam-dilution) of 
R(SiO) 65/54 = 0.66, since the source is unresolved by 
the SMT beams of 34" and 2874 in these lines. Using 
RADEX, we find that for Tkin = 15.5 K, the average H 2 
density in the SiO emission region is, Uav > 2 x 10® cm“® 


in order to fit the observed value of R(SiO)g 5 / 54 . The 
fractional SiO abundance is ~ 4x 10“®, close to the max¬ 
imum value assuming cosmic abundances (6 x 10“®), and 
comparable to the observed maximum circumstellar SiO 
abundance in oxygen-rich AGB stars {few x 10“®: Sahai 
& Bieging 1993). If the SiO emission comes from shocked 
gas at a higher kinetic temperature than that for CO, the 
minimum density required is higher. 

4.2. A Central Accretion Disk? 

Evidence of very fast outflows in 108005 comes from 
three independent probes: (i) CO data (presented here), 
(ii) Ha spectrum in SCetal08, and (iii) OH maser emis¬ 
sion (VLA mapping) reported in Zijlstra et al. (2001). 
The OH maser features cover the velocity range of 0 < 
Visr (km s“^) < 100, i.e., roughly the range covered by 
red half of the CO J=2-l profile. In contrast to OH 
and CO emission, the Ha absorption probes atomic gas. 
Since the Ha absorption feature is not spatially resolved 
in the ground-based long-slit spectra, it is difficult to di¬ 
rectly establish its relationship to the fast molecular out¬ 
flow, however the close agreement between the terminal 
outflow velocities derived for these indicates that they 
are closely associated. The atomic gas may be located 
inside the lobes and may constitute unshocked material 
of an underlying jet, and/or the interface between the lat¬ 
ter and the lobe walls, as proposed for IRAS 22036±5306 
(Scetal06). The jet-like feature seen in the HST image 
(Fig.[Tl inset) may represent the precessing jet’s signature 
close to its launch-site. 

The high-speed bipolar outflow in 108005 may be 
driven by an accretion disk. Bakker et al. (1997) find nu¬ 
merous narrow, double-peaked, chromospheric emission 
lines from neutral and singly ionized metals, and propose 
that these might arise in an accretion disk. We conjec¬ 
ture that the broad Ha wings seen towards this object 
might arise as a result of Raman-scattering of Ly/3 emis¬ 
sion generated by such a disk - a process that produces 
an Ha profile with a width that is a factor 6.4 larger 
than the Ly/3 width and a A“® wing profile. The very 
wide {FWZI ~ 2400 km s“^) Ha line wings in 108005 
show this shape, furthermore a weak emission feature 
around 6830A with EWZI ~ 150km s“^, correspond¬ 
ing to Raman-scattering of the 1032A component of the 
O VI doublet at AA1032,1038 is also seen (SCetal08); the 
7088A feature that corresponds to the 1038A component 
is a factor 4 weaker, and too weak to be visible. 

Other line-broadening mechanisms include electron 
scattering and emission from a rotating disk. Arrieta & 
Torres-Peimbert (2003) find that electron scattering re¬ 
quires extreme densities (e.g., Ue > 10^® cm“® in M2-9) 
making that an implausible mechanism. Keplerian rota¬ 
tion in a disk around the central star is too low to directly 
account for the broad line-width in 108005 (and PPNe in 
general: Sahai et al. 2011b) - taking the radius for the 
central post-AGB star of 108005 to be i? ~ 50i?© (Sli- 
jkhuis et al. 1991) and a nominal stellar mass of IM©, we 
find the maximum rotation speed is Vmx < 56.5km s“^, 
too low for generating the extreme line-wings of the Ha 
profile, even with the factor 6.4 increase provided by 
Raman-scattering. We therefore conclude that the ac¬ 
cretion disk in 108005 must be around a much smaller 
star, e.g., a main-sequence companion. 
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The Raman-scattering likely occurs in the innermost 
regions of the fast neutral wind seen via its blue-shifted 
absorption feature signature in the Ha profile. In this 
case, following Sahai et al. (2011b), we scale from 
the minimum scattering column density, Ng = 
cm“^ needed to achieve the Raman conversion effi¬ 
ciency necessary for producing the very broad observed 
line-widths (see Fig. 1 of Lee & Hyung 2000), to find 
that dMg/dt > (0.35 - 3.5) x IO-^Mq yr-i(rt„/30 AU) 
(D/3kpc) (t42,p/100km s“^) (A^s/10^®cm“^), where 
is the radius of the neutral wind where the Raman- 
scattering occurs. Vexp is set equal to 0.5 x FWHM 
of the broad CO J=2-l line profile, and similar to 
the neutral outflow’s average outflow-velocity derived 
by SCetalOS. The mass-loss rate requirement is easily 
met in 108005 since OBS05 derive a mass-loss rate of 
4.4 X 1O“®M0 yr“^ for a collimated fast wind in the 24- 
780 AU region around the central star, 

5. DISCUSSION 

108005 is very similar to other extreme-outflow oxygen- 
rich PPNs, IRAS 22036-f5306 (Sahai et al. 2006) and 
HD 101584 (Olofsson et al. 2015, Sivarani et al. 1999) in 
the properties of its collimated fast outflow (expansion- 
velocity, scalar-momentum, and kinetic-energy), and its 
Ha emission profile (very broad wings, P-Cygni absorp¬ 
tion). But in striking contrast to these objects, it appears 
to lack mm and submm continuum emission from its cen¬ 
tral region - the ratio of its 1.3 mm to 60 micron flux is 
< 10“"^ (3cr), compared to 0.6 x 10“^ and 1.5 x 10“^ for 
HD 101584 and IRAS 22036-f5306, respectively. 

Both PPNs, as well as the “disk-prominent” sub-class 
of young post-AGB objects (dpAGB objects: Sahai et 
al. 2011c), which unlike PPNs, show little or no extended 
nebulosity, and have central stars that are radial-velocity 
binaries (e.g., van Winckel et al. 2008), emit relatively 
strong mm/submm continuum emission from their cen¬ 
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tral regions. This emission has been attributed to the 
presence of substantial masses of cool, mm-sized grains 
(Sahai et al. 2011c, Gielen et al. 2007, de Ruyter et al. 
2006). Their origin is not understood at present but is 
potentially a key probe of important mass-ejection pro¬ 
cesses occurring during the late-AGB and post-AGB evo¬ 
lutionary phases, especially those that lead to the forma¬ 
tion of large dusty equatorial disks or torii. For example, 
in HD 101584, where a binary companion has been found 
from radial-velocity variations (Bakker et al. 1997), Olof¬ 
sson et al. (2015) propose a scenarion in which both the 
formation of the central, equatorially-dense mass struc¬ 
ture and the collimated outflow result from a common- 
envelope (CE) event - but the latter does not release 
enough energy to drive the mass-ejection, and another 
mechanism augments or even dominates it. 

The apparent lack of a substantial mass of material 
in the equatorial waist of 108005 suggests that perhaps 
it did not undergo a GE event, and its collimated out¬ 
flow (and collimated outflows in PPNs generally) may be 
launched differently, possibly from an accretion disk as 
we have suggested earlier. 108005 is thus a key post-AGB 
object for further detailed study. For example, ALMA 
can be used to probe its compact central region, in order 
to search for weak mm/submm continuum emission that 
may be present but was below our sensitivity limit, and 
for the presence of gas (and its kinematics) associated 
with its dusty waist. 
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TABLE 1 

SMA Observational Parameters 


Lines 

Freq. Range 
[GHz] 

Aiy 

[kHz] 

Array (# Ants.) 

Time^ 

[hr] 

Beam, PA 
// ^ // 0 

[K] 

Epoch 

yy/mm/dd 

i2CO(2-1) 

continuum^ 

USB®; 228.85-232.85 

228.9-229.6,229.9-230.2 

230.8-232.9 

812.5 

compact (6) 

1 

4.99 X 2.73,20.8 

100-230 

13/12/11 

13C0,C1*0(2-1) 

SiO(5-4),v=0 

continuum® 

LSB®: 216.85-220.85 

216.9-217.0,217.3-219.4 







i2C0(3-2) 

continuum^ 

LSB®: 342.0-346.0 
342.0-345.4 

812.5 

extended (7) 

2 

0.85 X 0.56,33.5 

240-400 

14/04/13 

continuum^ 

USB®; 354.0-358.0 








^On-source integration time 

*^The frequency range for the full spectrometer passband covering line and continuum regions 
'^The frequency ranges used for extracting the continuum 


Zijlstra, A. A., Chapman, J. M., te Lintel Hekkert, P., Likkel, L., 
Comeron, F., Norris, R. P., Molster, F. J., & Cohen, R. J. 2001, 
MNRAS, 322, 280 
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Fig. 1.— HST color-composite (red: F675W, green: F555W, blue: F439W) image of 108005 (log stretch). Panel size is 3794 x 3794. 
The diagonal linear features are diffraction spikes due to the bright central source in the relatively long-exposure F555W image. Inset is a 
false-color view (linear stretch: colorbar shows counts/s from 0 to 210) of the central 174 x 174 region (dashed box) in F439W, highlighting 
the curved jet-like feature emanating from the central source. 
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Fig. 2.— Molecular-line emission from 108005. The SMT spectra (shown above = 0 axis) have been shifted vertically for clarity, 
after scaling as follows: SiO J=6-5 and J=5-4 (x4), J=4—3 and 2-1 (x2), & J=2-l (x 10). The velocity ranges of the HV, MV, 

Sz LV components are marked in the SMA panel. The SMA J=2—1 spectrum, shifted below Tmb = 0, shows the total flux divided 

by the SMT conversion factor of 44.5 Jy/K, then scaled x2, for comparison with the SMT spectrum. The apparent difference in the SMT 
and SMA CO J=2-l line-profiles at Visr ~ 100km s~^ is probably due to non-linear baseline structure that is not removed by the linear 
baseline subtraction. 
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Fig. 3.— SMA CO emission integrated over different velocity (Visr) ranges (a) J=2-l, {top) blue: -150,-85, red: 130,200; (middle) 
blue: -85,-30, red: 70,130; (bottom) blue: —30,25, red: 25,75. (b) J=3-2, blue: —90,20, red: 30,130. Contour levels for CO J=2-l 
(3-2) are: minimum=50 (46)%, step=10 (15)% of the peak intensities. The central source in the HST image (marked with a cross), is 
offset (0^041,0^4) relative to the phase-center (at 0,0 with J2000 coordinates RA=08:02:40.7, Dec=-24:04:43.0), and consistent with the 
bipolar-outflow center, within the combined 1-cr uncertainty, ~(±0(M,±0(^4), of the absolute HST and SMT astrometry. 











